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ABSTRACT 

Context. PG 1 159-035 is the prototype of the PG1 159 spectral class which consists of extremely hot hydrogen-deficient (pre-) white 
dwarfs. It is also the prototype of the GW Vir variables, which are non-radial g-mode pulsators. The study of PG1159 stars reveals 
insight into stellar evolution and nucleosynthesis during AGB and post-AGB phases. 

Aims. We perform a quantitative spectral analysis of PG 1 159-035 focusing on the abundance determination of trace elements. 
Methods. We have taken high-resolution ultraviolet spectra of PG 1 159-035 with the Hubble Space Telescope and the Far Ultraviolet 
Spectroscopic Explorer. They are analysed with non-LTE line blanketed model atmospheres. 

Results. We confirm the high effective temperature with high precision (T e s= 140 000+5000 K) and the surface gravity of log g =7. 
For the first time we assess the abundances of silicon, phosphorus, sulfur, and iron. Silicon is about solar. For phosphorus we find an 
upper limit of solar abundance. A surprisingly strong depletion of sulfur (2% solar) is discovered. Iron is not detected, suggesting an 
upper limit of 30% solar. This coincides with the Fe deficiency found in other PG1159 stars. We redetermine the nitrogen abundance 
and find it to be lower by one dex compared to previous analyses. 

Conclusions. The sulfur depletion is in contradiction with current models of AGB star intershell nucleosynthesis. The iron deficiency 
confirms similar results for other PG1 159 stars and is explained by the conversion of iron into heavier elements by n-capture in the s- 
processing environment of the precursor AGB star. However, the extent of the iron depletion is stronger than predicted by evolutionary 
models. The relatively low nitrogen abundance compared to other pulsating PG1 159 stars weakens the role of nitrogen as a distinctive 
feature of pulsators and non-pulsators in the GW Vir instability strip. 
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1. Introduction 

PG 1 159-035 is the prototype star of the PG1 159 spectral class, 
which is a group of 40 extremely hot (r eff =75 000-200000 K) 
hydrogen-deficient post-AGB stars. They represent the evolu- 
tionary transition phase between Wolf-Rayet type central stars 
of planetary nebulae and the hottest non-DA white dwarfs. Their 
H-deficiency is probably caused by a late He-shell flash, which 
has laid bare the intershell region of the precursor red giant (see, 
e.g., a recent review by Werner & Herwig 2006). The study of 
the chemical composition of PG1 159 stars is therefore of partic- 
ular interest, because it allows us to directly analyse the outcome 
of mixing processes and nucleosynthesis in AGB star interiors. 

PG 1159-035 (=GW Vir) is also the prototype of variable 
objects among the PG1159 stars, which also comprises a num- 
ber of [WC] central stars. They are non-radial g-mode pul- 
sators. Besides the Sun, PG 1 159-035 is probably the star most 
intensively studied by asteroseismic methods (e.g. Kawaler & 
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Bradley 1994). Exploration of the interior structure allows in- 
sight into previous AGB evolution stages and is thus of signifi- 
cance well beyond the PG1 159 stars themselves. The basic pul- 
sation driving mechanism is the k - y effect of C and O acting 
beneath the photosphere (Starrfield etal. 1983), however, many 
details concerning the role of abundances of other species are 
still debated (e.g., Gautschy etal. 2005). 

The particular significance of PG 1 159-035 for insight into 
(post-) AGB stellar structure and evolution motivates our work 
to perform a detailed spectroscopic analysis of this star. We 
have taken high-resolution UV spectra with the Far Ultraviolet 
Spectroscopic Explorer (FUSE) and the Hubble Space Telescope 
(HST) that cover the wavelength interval from the Lyman edge 
up to 1730 A. The emphasis of our analysis is on the identifica- 
tion and abundance analysis of trace elements. The results will 
be discussed in the framework of late He-shell flash and AGB 
evolution models as well as pulsation theory. 

PG 1 159-035 was discovered in the Palomar Green Survey 
(Green et al. 1986, Wesemael et al. 1985) and was first subject to 
a quantitative non-LTE spectral analysis by Werner et al. (1991). 
Based on optical spectra, r e ff=140000 K and log g =7.0 were 
derived and the main atmospheric constituents were found to be 
He=33%, C=50%, and 0=17% (mass fractions). Later work in- 
cluding the present one confirmed these parameters and we con- 
centrate here on other chemical elements. 

UV spectroscopy has been performed with the International 
Ultraviolet Explorer (IUE, Liebert etal. 1989), the Hopkins 
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Table 1. Identified photospheric lines (rest wavelengths) in the 
FUSE spectrum of PG 1 159-035. 
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Ultraviolet Telescope (HUT, Kruk & Werner 1998), as well as 
with the Faint Object Spectrograph aboard HST (Werner & 
Heber 1993, Dreizler & Heber 1998). It revealed that, as in the 
optical band, the spectra are dominated by lines from He n, C iv 
and O vi. The resolution and S/N of the spectra presented here 
are far superior. We will refer to results from earlier observations 
where relevant. 

The paper is organized as follows. We present our new 
HST and FUSE observations in Sect.|2]and identify photospheric 
lines. In Sect.|3we describe our model atmosphere and line for- 
mation calculations. Sect. 0] contains the results of the spectral 
analysis and we conclude with a discussion of our results in 
SectJU 



Table 2. Photospheric lines (rest wavelengths) in the HST spec- 
trum. 
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Fig. 1. The complete FUSE spectrum of PG 1159-035 compared to our final model spectrum (black graph). The detected pho- 
tospheric lines are identified. Labels in brackets denote lines that are visible in the model but not unambiguously seen in the 
observation. Overplotted is the same model but now attenuated by ISM absorption lines (grey graph; in the online version of this 
paper the graph is red). The ISM lines are not labeled. Photospheric model parameters are given in Table Q In all figures of this 
paper the observed spectra are shifted such that the photospheric lines are located at rest wavelengths. 



2. Observations 

2. 1 . FUSE and HST spectroscopy 

The FUSE spectrum (observation id Q1090101) was obtained on 
May 10, 2001, using the LWRS aperture. The observation was 
split evenly into three exposures in successive orbits for a total 
exposure time of 6324 sec. The spectrum covers the range from 
the Lyman edge up to 1 187 A, i.e., there is a narrow range that 



overlaps with the HST data. The data were reduced with ver- 
sion 3.1.7 of the CalFUSE pipeline. The spectral resolution was 
rather poor in this first reduction, about 10000, corresponding 
to «0. 1 A, indicating that motion of the optics during the expo- 
sure might have been significant (Sahnow etal. 2000). The ex- 
posures were then reprocessed in segments of 400-500 sec each. 
The spectra from these separate sub-exposures were coaligned 
by cross-correlation on narrow interstellar absorption lines, and 
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Fig. 2. Similar to Fig.^ Here we show the complete HST spectrum of PG 1 159-035. In addition to the photospheric lines we also 
label the interstellar components of the resonance doublets from C iv and Si iv. 
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Table 3. Radial-velocity shift of selected photospheric lines. 



Ion 


/W / A 


■^observed / A 


v rad / km/s 


Nv 


1238.82 


1239.06 


58.08 


Nv 


1242.80 


1243.04 


57.89 


Sirv 


1393.76 


1394.07 


66.68 


Civ 


1548.19 


1548.48 


56.16 


Civ 


1550.77 


1551.06 


56.06 


Ov 


1371.30 


1371.57 


59.90 



0.6 



Mean 



59.13 ±1.62 



combined. The resulting spectral resolution was much closer to 
the expected value of 20 000. 

The spectra from the four channels were then combined to 
produce a single spectrum. The spectra were first resampled onto 
a common wavelength scale: the pixel width of 0.013 A was not 
changed, but non-integer pixel shifts were required. This causes 
a slight smoothing of the spectra, but the effect is small in com- 
parison with the instrumental resolution width. At some wave- 
lengths additional shifts of 1 to 2 pixels were required to coalign 
the spectra from different channels, as a result of small resid- 
ual distortions in the detector pixel scale. Regions of known bad 
pixels were discarded. The broad depression in the LiF lb spec- 
trum induced by the grid wires in the detector was corrected by 
normalizing the LiF lb spectrum to LiF 2a with a 100-pixel slid- 
ing boxcar filter. A similar but weaker depression in the SiC lb 
flux over wavelengths of 977-992 A was corrected by normal- 
izing to the SiC 2a spectrum in the same manner. The spectra 
were weighted by the statistical uncertainties on a pixel-by-pixel 
basis when combined to produce the final spectrum, which is 
displayed in Fig.^ 

The HST spectrum was obtained on June 9, 2004, with the 
Space Telescope Imaging Spectrograph (STIS) using the grat- 
ing E140M and the 0.2" x 0.2" aperture. The exposure time 
was 2275 sec, divided evenly among three exposures. The spec- 
trum covers the range 1 150-1730 A with a resolution of 90 000 
corresponding to roughly 0.02 A. The standard calibrated data 
files provided by the Multimission Archive at Space Telescope 
(MAST) were used, which were reduced using the most recent 
version of the CalSTIS pipeline (v 2.20) and associated reference 
files. These spectra were resampled onto a common wavelength 
scale and combined, weighting by the statistical uncertainties of 
the data. The pixel width of the wavelength scale varied linearly 
from 0.0123 A at 1140 A to 0.01826 A at 1729 A, closely ap- 
proximating the sampling of the individual echelle orders. Pixels 
at the ends of the echelle orders with bad quality flags were dis- 
carded. The final merged spectrum is displayed in Fig.|5] 

The FUSE and STIS spectra in Fig s.^ and |3 are overplotted 
with the astrophysical flux (cgsA units) from our final model, 
normalized to the local continuum. 



2.2. Line identification 

Line identification was guided by model-atmosphere spectra that 
involve atomic data from different sources, which are described 
in detail in Sect. 13.11 

Table ^ lists al1 photospheric lines detected in the FUSE 
spectrum. The majority stems from He it, C rv, and O vi and was 
already detected in the HUT spectrum (Kruk & Werner 1998). 
However, thanks to its higher resolution the FUSE spectrum al- 
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Fig. 3. The O v 1371.30 A line compared with two models with 
different effective temperatures. 



lows for the identification of individual multiplet components of 
several C iv and O vi lines. The discovery of highly ionized neon 
and fluorine lines (Ne vn 973.33 A, F vi 1 139.50 A) in the FUSE 
spectra of PG 1 159-035 and other PG1 159 stars was announced 
recently (Werner etal. 2004, 2005). In addition, we present here 
the discovery of the S vi 933.38/944.52 A resonance doublet. We 
also identify for the first time several lines from Si v; we discuss 
this further below together with other silicon features (Si iv-vi) 
seen in the HST spectrum. 

We cannot confirm our earlier claimed detection of N v lines 
in the HUT spectrum. A broad and shallow feature at 1050 A 
was attributed to the 4d-6f and 4f-6g transitions at 1048.3 and 
1049.7 A, which was supported by a model spectrum calculated 
with 1% nitrogen. In this paper we derive from the Nv reso- 
nance doublet in the HST spectrum that the N abundance is a 
factor of 10 lower. As a consequence, the two lines at 1048.3 and 
1049.7 A are very weak in our final model (Fig.[Q. The two nar- 
row lines in the FUSE spectrum close to the wavelengths of the 
N v lines are interstellar lines of Ar 1 1048.22 A and H 2 R(0)4_o 
1049.374 A. As in the HUT spectrum, the FUSE spectrum also 
shows a broad feature at 1048-1052 A. An alternative identifi- 
cation could be a very highly excited C iv line (n=4-l 1) which, 
however, cannot be checked with our models because the upper 
level is not treated in NLTE in our model atom. 

The FUSE spectra show that the C 111 977.0 A line detected in 
the HUT spectrum is entirely of interstellar origin: the Doppler 
shift matches the velocity of the other interstellar absorption fea- 
tures, and the predicted line feature in our model atmosphere 
spectrum is very weak. 

We will later discuss the non-detection of the Pv 
1117.98/1128.01 A resonance doublet. It was discovered in 
cooler PG1 159 stars (Reiff et al. 2006) but PG 1 159-035 is ob- 
viously too hot to exhibit these lines, so that we will only be able 
to derive an upper abundance limit. 

Table|5]lists all photospheric lines identified in the HST spec- 
trum. The majority stems from C iv and O vi. In addition we see 
a broad He 11 1640 A line, blended in a wide trough with C iv and 
O vi lines (Fig.[2J- The O v 1371 A line is seen with a narrow and 
deep profile (Fig.[3}. Concerning the N v 1239/1243 A resonance 
doublet, our spectrum allows for a clear distinction between in- 
terstellar and photospheric components and, thus, for a reliable 
abundance determination. 
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We also discovered the Si iv 1394/1403 A resonance doublet, 
but it is very weak (Fig.[6}. Identification of the blue component 
is rather doubtful although it should be stronger than the red one. 
We were able to find stronger silicon lines from higher ionisation 
stages, Si v and Si vi, in the FUSE and HST spectra (Fig.0. To 
our knowledge, this is the first detection of such lines in a stellar 
photospheric spectrum. Accordingly, Fig.0]shows that these two 
ionisation stages are more strongly populated than Si iv through- 
out the atmosphere. 

There is a quite strong absorption feature at 1319.7 A (low- 
est panel of Fig.0. It was discovered by Feibelman (1995) in 
an IUE spectrum of the PG1159-type central star of NGC246. 
A close inspection of the profile shows that the feature is a 
blend of two lines, separated by about 0.2 A. We suggest that 
they are lines from neon and silicon although our model is nei- 
ther able to fit the position accurately nor to match the strength 
of the feature. One of the two lines could be a component 
of the Nevn 2p 1 P°-2p 23 P intercombination triplet located at 
1319.78/1335.35/1344.52 A. (Unfortunately, the two other com- 
ponents are blended by a strong interstellar line and a Civ 
line, preventing an unambiguous identification.) The second line 
blending with this Ne vn component is probably Si v 3s 'P - 
3p'D, located at 1319.60 A. 

We do not confirm the claimed detection of an Ovn line 
(1522 A) in the IUE spectrum (Feibelman 1999). 

Probably most other absorption lines in the HST spectrum 
are of interstellar origin, however, there is a number of features 
that we think are photospheric but which we are unable to iden- 
tify (see Sect.l2~51. 

From photospheric lines in the HST spectrum we have de- 
termined a radial velocity shift of v ra ^ = 59.13 + 1.62 km/s 
(Tab. Using the same set of lines in IUE spectra, Holberg 
etal. (1998) determined v ra </ = 50.10 +1.10 km/s. For the same 
set of ISM lines that was used by Holberg etal. (1998) we de- 
rive v ra d = -9.78 + 0.59 km/s as compared to their value of 
v ra rf = -17.73 + 1.06 km/s. In all our Figures showing the HST 
spectrum we have shifted the observation such that the photo- 
spheric lines are at their rest wavelength position. 

All observed spectra displayed here are convolved with 
Gaussians with FWHM=0.05A and the synthetic spectra with 
FWHM=0.1A. 

2.3. Unidentified photospheric lines 

In Table |4] we list unidentified lines that we think are photo- 
spheric. We have carefully compared the FUSE and HST spectra 
of PG 1159-035 to spectra of other hot post-AGB stars whose 
photospheric and ISM radial velocity systems have different sep- 
arations. In the Table we note in which objects we have seen the 
same feature as in PG 1159-035. Most useful for this compar- 
ison were other HST/STIS spectra, namely that of the PG1159 
star NGC 7094, which we observed with the same instrumental 
setup as PG 1159-035, and a very high resolution spectrum of 
NGC 246 (Jenkins priv. comm.; Rauch etal. in prep.). We also 
compared the PG 1 159-035 spectra with FUSE and HST spec- 
tra of various hydrogen-rich central stars, which we are currently 
analysing (Traulsen et al. 2005, Hoffmann et al. 2005). 

The unidentified lines likely stem from metals other than 
CNO, whose line lists we consider complete. They are probably 
not unknown iron lines, because we do not detect the strongest 
known lines of this species (Sect. l4.lt . We suggest that the lines 
stem from highly ionized light metals like neon or magnesium 
for which line lists are inaccurate and incomplete. Some of the 



Table 5. Summary of the model atoms used in the model atmo- 
sphere and line formation calculations. The numbers in brackets 
give the individual line numbers summed into superlines for the 
heavy metal ions. 
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IX 
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25 


(36 843) 




X 
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unidentified lines are quite strong and the most prominent one is 
located at 1270.25 A. 

3. Model atmospheres and synthetic spectra 

We used the Tubingen Model Atmosphere Package tmap 
(Werner et al. 2003, Rauch & Deetjen 2003) to compute plane- 
parallel non-LTE line blanketed model atmospheres in radiative 
and hydrostatic equilibrium. The models include the most abun- 
dant elements: He, C, N, O, Ne. At first, we computed a small 
grid varying the elemental abundances as well as r ei f and log g . 
The comparison of the resulting line profiles to the HST and 
FUSE spectra confirmed previously published values (except for 
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Table 4. Unidentified photospheric lines in the FUSE and HST spectra (rest wavelengths). We remark whether the lines are de- 
tectable in high-resolution spectra of other PG1 159 stars or H-rich central stars as well. 

wavelength / A Remark Possible identification 

927.30 
939.57 
965.35 
996.58 
999.42 
1000.08 
1006.20 

1041-1047 some features 

1052.00 
1055.66 
1057.74 

1 1 1 0- 1 1 20 some features 

1125.68 
1127.57 
1129.5 
1131.17 
1131.80 
1139.04 
1141.60 
1148.42 

1148-1152 Kl-16, RXJ2117+3412, Lo4, PG1144+005, PG1707+427 broad dip, Civ 4-11 

1 155-1 160 some features 

1264.15, 1264.95, 1265.95 NGC 7094, H-rich CSPN 

1267.65 NGC 7094, H-rich CSPN 1267.64 Si v 3s»P°-3p 3 P, no f-value available 

1270.25 NGC 246, NGC 7094, H-rich CSPN 

1298.95 H-rich CSPN 

1307.60 H-rich CSPN 

1356.00 H-rich CSPN 

1366.40 RX J2 117.1 +3412 (HST/GHRS, Werner etal. 1996) 

1406.45 H-rich CSPN 

1480.20 NGC 246, not in H-rich CSPN 

1486.75 H-rich CSPN 

1645.00, 1645.60 NGC 7094, H-rich CSPN 

1679.55 NGC 7094, H-rich CSPN 



the N abundance). In particular, the Ov 1371 A line is a very 
sensitive r e ff indicator and we confirm our result from an earlier 
fit to this line taken with HST/FOS (Werner & Heber 1993), but 
now with a high precision: 140000+5000 K. In a second step 
we adopted these parameters for T e g, log g , and abundances 
of He, C, O, Ne and included several other elements to deter- 
mine their abundances. For this, we performed line formation 
iterations keeping fixed the atmospheric structure. In this man- 
ner we investigated the line profiles of N, Si, P, S, and Fe. The F 
abundance was determined recently (Werner et al. 2005) and this 
value is adopted for our final spectrum synthesis calculation. 



3.1. Atomic data 

Table [5] summarizes our model atoms. We used our standard 
model atoms for He, Ne, F, Fe that were described in previously 
cited work. Other model atoms are extended versions of previous 
datasets or designed completely a new. We briefly summarize 
these advancements. We used the databases of NIST 1 (National 
Institute of Standards and Technology) and CHIANTI 2 for 



1 http://physics.nist.gov/ PhysRefData/ASD/index.html 

2 http://wwwsolar.nrl.navy.mil/chianti_direct_data.html 



the energy levels and the Opacity Project and Iron Project 
(TIPTOPbase 3 ) for the oscillator strengths. 

The model ions of C iv, N v and O vi were updated and com- 
pleted. The most important improvement was made for O vi. The 
FUSE and HST spectra of PG 1 159-035 and of other hot stars 
reveal that the positions of many listed O vi lines are inaccurate. 
We are therefore able to improve atomic data by determining the 
line positions from observations. In Table|S]we list for all these 
lines the literature values and our new observationally derived 
positions. We find deviations up to 0.6 A. In our spectrum syn- 
thesis we apply the new line positions. 

The model atom for phosphorus is entirely new while for sul- 
fur we use a model designed by Miksa et al. (2002). The silicon 
model atom was extended to include lines from ionisation stages 
higher than Si iv. 

Figure@]displays the ionisation stratification of all elements 
considered in the final atmospheric model. 

4. Spectral analysis 

4.1. Abundance determination of N, Si, P, S, and Fe 

Nitrogen. — Our HST spectrum allows us to distinguish the inter- 
stellar and photospheric components of the N v resonance dou- 

3 http://vizier.u-strasbg.fr/topbase/ 
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Fig. 4. Depth dependence of temperature, electron density, and ionisation fractions of all considered chemical elements in the final 
atmospheric model (for its parameters see Tab.0. 
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Table 6. Newly determined positions of O vi lines as used in our 
improved model atom. We compare them to the literature values. 
The measurements were made in the PG 1159-035 spectra de- 
scribed in this paper (J06) as well as in spectra of other PG1 159 
stars (Reiff et al. 2006, R06) and hot central stars (Traulsen et al. 
2005, T05). 



1.0 



Transition 



/ioid / A 



, / A Remark 



4p 2 P? /2 
4p 2 P° /2 

4p 2 P« p - 
4p 2 P° /2 - 
4p 2 P° /2 - 

4f 2 F° 5p 
4f 2F S/2 



4f 2 F 



7/2 



5p 2 P° p 
5p 2 P° p 
5p 2 P° 



-V2 



5s 2 S I/2 
-5s 2 S 1/2 

■ 5d 2 D 3/2 

■ 5d 2 D 3/2 

■ 5d 2 D 5/2 

5d 2 D 3/2 
5d 2 D 5/2 
5d 2 D 5/2 

■ Vd 2 D 3/2 

■ 7d 2 D 3/2 

■ 7d 2 D 5/2 



5p 2 P° p 
5p 2 P° p 
5p 2 P° 



3/2 



-V2 



8d 2 D 

8d 2 D 3/2 
8d 2 D 5/2 



5/2 



5f 2 F° 
5f 2 F" 



5/2 



5d 2 D 3/2 - 7f 2 F 
5d 2 D 5/2 - 7f 2 RJ 
5d 2 D 5/2 - 7f 2 F^ 2 



7d 2 D 3/2 
7d 2 D 5/2 
7d 2 D 5/2 

■ 7g 2 G 7/2 

■ 7g 2 G 7 / 2 
7g 2 G 9/2 

-7f 2 F? 



5f 2 F 



5f 2 F° 5/ 2 - 
5f 2 F° /2 - 

5f 2 F° - 

Ji r lp 

5g 2 G 7/2 
5g 2 G 7/2 
5g 2 G 9/2 

5g 2 G 7/2 
5g 2 G 9/2 
5g 2 G 9/2 • 



7f 2 F° 



5/2 



7f 2 F° 



7/2 



7/2 



• 7h 2 H° 



■ 7h 2 H° 



9/2 



■ 7h 2 H° 



9/2 



1 1/2 



1171.56 
1 172.44 

1080.88 
1081.62 
1081.52 

1126.46 
1126.35 
1126.46 

1262.10 
1262.20 
1261.82 

1018.21 
1018.28 
1018.15 

1289.82 
1289.97 
1289.94 

1293.14 
1292.74 
1292.77 

1291.91 
1291.94 
1291.91 

1292.12 
1292.09 
1292.12 

1291.99 
1292.02 
1292.02 



1171.10 
1172.00 

1080.60 
1081.34 
1081.24 

1126.28 
1126.17 
1126.28 

1261.70 
1261.80 
1261.42 

1018.81 
1018.88 
1018.75 

1290.06 
1290.21 
1290.18 

1293.04 
1292.64 
1292.67 

1291.81 
1291.84 
1291.81 

1292.02 
1291.99 
1292.02 

1291.89 
1291.92 
1291.92 



T05 
T05 

R06 
R06 
R06 

R06 
R06 
R06 

J06 
J06 
J06 

J06 
J06 
J06 
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T05 
T05 

T05 
T05 
T05 

T05 
T05 
T05 

T05 
T05 
T05 

T05 
T05 
T05 



blet (Fig. |3- We find a good fit with N=0.001 (mass fraction). 
This is one dex lower than the value derived by Dreizler & Heber 
(1998) from the HST/FOS spectrum that could not separate the 
interstellar and photospheric components. 

Silicon. — Figures [5] and Q show the resonance doublet of 
Sirv as well as lines from Siv-vi together with two synthetic 
spectra with solar and half-solar Si abundance. We conclude that 
the Si abundance is about 0.5 solar. 

Phosphorus. — Figure [8] shows the FUSE spectral region 
where the resonance doublet of Pv is located. The doublet is 
not detected. By comparison with models having solar and ten 
times solar P abundance we conclude that the P abundance is at 
most solar. 

Sulfur. — Figure [5] shows the resonance doublet of S vi to- 
gether with two model calculations with solar and 2% solar 
abundance. It is obvious that the S abundance is strongly de- 
pleted. We find a best fit at 2% solar. 




1238 



1240 1242 
wavelength / A 



1244 



Fig. 5. The N v resonance doublet compared with two models 
with different nitrogen abundance. The model with N=0.001 
(mass fraction) fits the observed profiles very well. Note the 
weak blueshifted ISM components which are included in the 
model. 



Iron. — We have closely inspected the HST and FUSE spec- 
tra looking for iron lines. The only Fe lines predicted from our 
models stem from Fevn. The dominant ionisation stages are 
Fevm-ix (Fig.@} but no detectable lines with accurately known 
wavelength position are known from this ion (i.e. no such lines 
are listed in the respective Kurucz (1991) POS files in the rel- 
evant wavelength region). The Fe vn lines are too weak and we 
are not able to make a clear identification. From a comparison 
of the PG 1 159-035 spectrum with those of (cooler) H-rich cen- 
tral stars and with our model we found that the Fevn line at 
1332.38 A should be among the strongest. There is an absorp- 
tion feature in the HST spectrum at this location CFig. 1101 . but 
the case is less convincing with other lines. We cannot detect be- 
yond doubt any iron line. The computed line profiles from mod- 
els with solar and 0. 1 solar Fe abundance suggest an upper limit 
of Fe^0.3 solar. 

We have also tried to look for lines from nickel, which is an 
important element in the context of the Fe depletion discovered 
in PG 1159-035 as well as in other PG1159 stars. The situa- 
tion for Ni is even less certain than for Fe. The highest ionisa- 
tion stage having precisely known line positions in the HST and 
FUSE ranges is Ni vi. The ionisation structure of Ni in the atmo- 
sphere is comparable to Fe and our models predict no detectable 
Ni vi lines at all because of the very high effective temperature. 

4.2. Absorption by interstellar material 

Absorption by interstellar gas was modelled to facilitate analysis 
of the photospheric spectrum. Species detected in the data and 
included in the model are: H i, D i, C i, C n, C n*, C m, C iv, N i, 
N n, N m, N v, O i, Mg n, Al n, Si n, Si m, Si iv, S n, S m, S iv, Ar i, 
Mnn, Fen, and H2 (J=0,l,2,3). All interstellar features in the 
FUSE and STIS spectra were fitted simultaneously. A satisfac- 
tory fit to all the neutral and singly-ionized species typical of the 
warm diffuse ISM was obtained with a single component at a he- 
liocentric velocity of -9.7 km/s. The C 1 and H2 features were fit 
with a separate component; the Doppler width was much smaller 
than that of the warm gas, as expected, but the radial velocity was 
the same. The C m, Si m, S m and S iv absorption features were 
also at about the same velocity as the low-ionization species, 
with some scatter of a few km/s. The highly-ionized species of 
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Fig. 6. The photospheric Siiv resonance doublet is barely de- 
tectable while the ISM component (included in the model) is 
relatively strong. The two model profiles suggest that the silicon 
abundance is about 0.5 solar. 
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Fig. 7. Examples for newly discovered lines from highly ionised 
silicon (Si v-vi) compared with models with solar and 0.5 solar 
Si abundance. 
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Fig. 8. The non-detection of the P v resonance doublet and the 
two model spectra suggest that the P abundance is at most solar. 
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Fig. 9. The relative weakness of the S vi resonance doublet sug 
gests a strongly subsolar sulfur abundance. 



Si iv, C iv, N v, and O vi exhibit a more extended distribution in 
velocity. The N v profiles show two distinct components, one at 
-12 km/s and one at +14 km/s. The Si iv, C iv, and O vi profiles 
span the same range of velocity, but the substructure is not as 
clearly resolved. These velocities range from a value similar to 
that of the ISM absorbers to one that is intermediate between the 
ISM material and the photosphere. This suggests that the highly- 
ionized material is circumstellar; a detailed analysis of this ma- 
terial is outside the scope of this paper. The final ISM model 
was evaluated for each pixel in the observed spectra, and multi- 
plied by the photospheric model to produce the plots shown in 
Figures^and|2 

We derived the neutral H column density towards 
PG 1159-035 by a fit to the interstellar Ly ff line (Fig.0. Prior 
to fitting, the spectrum was normalized by the photospheric 
model to account for the broad Hen and Civ absorption fea- 
tures; consideration of the photospheric Civ lines at 1210.61 A 
(4s-7p) and 1213.56 A (4d-8f) is particularly important when 
attempting to derive «h from Ly ff for stars in this spectral class. 
Uncertainty in the proper background level was assessed by per- 
forming a series of fits to the Ly„ profile in which the residual 
background level was varied from -0.06 to +0.04, in steps of 
0.02. The best fit was obtained with a background level of 0.0 
and«n = 1.50T0 20 crrT 2 . The mean background level in the cen- 
tral portion of the Ly,, trough, apart from the narrow geocoronal 
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Fig. 10. Computed Fe vn lines (with solar and 0. 1 solar Fe abun- 
dance) compared with observations. There is no clear evidence 
for the presence of iron lines in PG 1 159-035. This suggests that 
iron is depleted and we estimate an upper abundance limit of 0.3 
solar. 



emission feature, is w-0.02. For this value of the background, the 
best-fit n H is 1.45-10 20 cm -2 ; we therefore take 0.05-10 20 as rep- 
resentative of the uncertainty in the determination of «h arising 
from the uncertainty in the background determination. 

The effect of the uncertainty in T e ff was assessed by repeating 
the series of fits on spectra that had been normalized by models 
with r eff =136000K and 146 000 K. The resulting chi-squared 
contours had the same minimum position and same gradients 
away from the minimum, so uncertainty in T e g does not appear 
to contribute significantly to the overall uncertainty in «h- 

Our adopted value of njj = 1.50 + 0.05 • 10 20 cm -2 is fairly 
consistent with earlier measurements («h = 1-0- 10 2() cm -2 , Kruk 
& Werner 1998). We found that it is not necessary to apply inter- 
stellar reddening to our model in order to reproduce the contin- 



uum shape of the observed spectrum. This is also in accordance 
with what we found from the HUT data (Kruk & Werner 1998). 

4.3. Summary of spectral analysis 

TableQsummarizes the results of our spectral analysis including 
results from earlier works. We also list mass, luminosity, and 
distance estimates taken from Werner et al. (1991). 

The overall fit of our final model to the FUSE and HST 
spectra is good (Figs.^and[2j. Let us comment on those few 
photospheric features where deviations between model and ob- 
servations are obvious. The Civ 4p-10d line at 1060 A is not 
seen in the observed spectrum and is too strong in the model. 
The explanation could be that the upper level population is not 
computed realistically enough because the n — 10 levels are the 
highest ones treated in NLTE. The n > 10 levels are treated 
in LTE which might be insufficient. The computed profiles of 
two Ovi lines have no detectable counterpart in the observed 
spectrum, namely the 5p-7s and 5d-7p transitions at 1303.8 and 
1302.5 A. As we have shown above, the exact wavelength po- 
sition of many O vi transitions, all of them involving the n — 5 
levels (see Table [6j, are uncertain within a few tenths of an A. 
For the two lines in question we could not find plausibly shifted 
counterparts in the observed spectrum. It could be that they are 
blended with the close-by strong interstellar features. 

Concerning the continuum shape, our model overestimates 
the flux in the long-wavelength region of the HST spectrum by 
a few percent when normalized to the short-wavelength region. 
This means that the HST spectrum is harder than our model, 
which is unphysical. Because of the high r e fj- our model has a 
Rayleigh-Jeans shape and the stellar spectrum cannot be steeper. 
Hence, there are problems with the flux calibration of the HST 
spectrum. On the compressed scale of Fig.|5J the discrepancy is 
only apparent at the longest wavelengths where it has grown to 
nearly 10 percent, but on an expanded scale the effect appears to 
begin longward of about 1420 A. 

5. Summary and discussion 

We have analysed high-resolution FUSE and HST spectra of 
PG 1 159-035 covering the wavelength range 912-1730 A. They 
are dominated by lines from highly ionized carbon and oxygen 
plus several light metals (Tables [2 and |2ji. More than forty rela- 
tively strong absorption lines in the UV spectra remain uniden- 
tified (Table@J. They most probably stem from metals between 
the CNO elements and the iron group. We discovered that many 
spectral lines from the O vi ion have inaccurately known wave- 
length positions (up to 0.6 A). Our high resolution spectra allow 
for a precise wavelength determination (< 0.05 A) for about 30 
lines (Table|6ji. 

The principal aim of this work was the abundance analysis 
of trace elements beyond the main atmospheric constituents He, 
C, and O (Table0. Since the PG1 159 stars display former inter- 
shell matter on their surface as a consequence of a (late) He-shell 
flash, the elemental abundances allow us to place constraints on 
nucleosynthesis and mixing processes taking place in the AGB 
phase. We summarize and discuss our results in brief. For a more 
general discussion of abundances in PG1 159 stars in this context 
see Werner & Herwig (2006). 

Nitrogen. — For the first time we succeeded resolving the 
photospheric and interstellar components of the N v resonance 
line, allowing a reliable abundance determination. We find 
N=0.001 (mass fraction), which is an order of magnitude 
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Table 7. Photospheric and stellar parameters of PG 1159-035. 
We list the direct results from spectroscopic analyses and de- 
rived quantities (mass, luminosity, distance) from comparison 
with evolutionary tracks. Elemental abundances are given in 
mass fractions (2nd column) and relative to solar abundances 
(Grevesse & Sauval 2001 values; 3rd column). 



Parameter 


Result 


Abundance s 
(solar units) 


Ref. 


7k/ K 


140000 ±5000 




(1), this work 


log g (cgs) 


7.0 ±0.5 




(1) 


H 


< 0.02 


< 0.027 


(2) 


riC 


U.J J 


1 ^ 


v- ) 


l_ 




1 70 




N 


0.001 


1.2 


this work 


O 


0.17 


22 


(3) 


F 


3.2 • 10~ 6 


6.3 


(3) 


Nc 


0.02 


11 


(3) 


Si 


3.6 ■ 10- 4 


0.5 


this work 


P 


< 6.4 • lO" 6 


< 1 


this work 


s 


5.0 • 10~ 6 


0.02 


this work 


Fc 


< 3.9 • 10~ 4 


< 0.3 


this work 


«h / crrr 2 


1.5- 10 20 




this work 


E(B-V) 


0.0 




this work 


M/M Q 


0.60!;" 




(1) 


L/L 


2J3+ 'I 




(1) 


d/kpc 


0.80^ 




(1) 



References: 
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less than the upper limit previously known. We conclude that 
PG 1159-035 suffered a late thermal pulse (LTP), i.e., a late 
He-shell flash during the previous H-burning post-AGB phase. 
The flash did not occur on the white dwarf cooling track (a very 
late TP, VLTP), because in this case a larger N abundance is ex- 
pected. The N abundance in PG 1 159-035 is about one tenth of 
that found in other GW Vir pulsators. Hence, the presence of a 
high N abundance does not strictly correlate with the occurrence 
of pulsations as previously found (Dreizler & Heber 1998). 

Fluorine. — The strong F overabundance in PG 1 159-035 as 
well as other PG1 159 stars was found and discussed in previous 
work (Werner et al. 2005). It emphasizes the important role of 
AGB stars as effective Galactic fluorine producers. 

Neon. — The Ne overabundance in PG 1159-035 and other 
PG1 159 stars was also discussed previously (Werner et al. 2004). 
It confirms the strong convective overshoot occurring during He- 
shell burning in thermally pulsing AGB stars. 

Silicon. — We were able to determine the Si abundance in 
PG 1159-035 for the first time. This is based on the Siiv res- 
onance doublet and on several lines from Si v-vi. The latter have 
never been observed before in any stellar photospheric spectra. 
Within the error limit Si is solar. This is in accordance with ex- 
pectations from stellar evolution models. 

Phosphorus. — We found a solar abundance as an upper 
limit. Current evolutionary models, however, predict overabun- 
dances of 4-25 times solar. 

Sulfur. — Based on our first discovery of the S vi resonance 
doublet we find an extreme underabundance of S=0.02 solar. 
This strongly contradicts evolutionary models that predict only 
slight depletions (0.6 solar). One might speculate that the S de- 
pletion found here is related to the sulfur abundance anomaly in 
Type II (i.e. C-rich) planetary nebulae. Many of these exhibit an 



unexpectedly low S/O ratio in contrast to many Type I PNe and 
H II regions (Henry et al. 2004, 2006). It was argued, however, 
that this anomaly in PNe is not real but due to the use of a model- 
dependent (and for some unknown reason incorrect) ionisation 
correction factor. 

Iron. — We determine an upper limit of Fe=0.3 solar. This is 
in agreement with results for three other PG1 159 stars where we 
found upper limits of 0.01-0. 1 solar (Miksa et al. 2002). We have 
suggested that Fe was transformed into Ni and heavier elements 
by s-process neutron captures; however, such strong Fe depletion 
in the intershell is not predicted by current AGB star models. On 
the other hand, this idea is confirmed by the strongly subsolar 
Fe/Ni ratio found in Sakurai's VLTP star (Asplund et al. 1999). 
In the case of PG 1159-035 we are not able to determine any 
constraint for its Ni abundance. 

While the abundances of the main atmospheric constituents 
of PG1 159 stars can be explained quantitatively by evolutionary 
models for AGB stars that later go through a (late) thermal pulse, 
a close look at trace elements abundances reveals that many de- 
tails remain unclear. The extreme depletions of sulfur and iron 
in PG 1 159-035 are most unexpected. 
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